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ABSTRACT 


J o 

Recent measurements of the electron temperature (T % (1.5±.5) x 1CT ~K) 

6 

2 

show that 3 * nk(Te + Tp)/(B /8n)*l, indicating that magnetic and thermal 
pressures are equally important in the solar wind. Hydromagnetic theory, 
rather than hydrodynamics, must thus be used. MHD theory predicts sr /eral 
types of discontinuities. Tangential and rotational discontinuities are 
important in cosmic ray scattering ther y. Observations show that most 
discontinuities in the interplanetary magnetic field are tangential. For- 
ward and reverse fast shocks and slow shocks, predicted by hydromagnetic 
theory, have all been observed. Forward fast shocks can accelerate inter- 
planetary particles, increasing their energy several fold to ~ 1 MeV. 

Fermi acceleration of particles trapped between the shock and magnetic 
fluctuations can explain the observations. 






I. Introduction 


The aim of this paper is to discuss recent work by the author and 
his colleagues at Goddard Space Flight Center concerning discontinuities 
in the interplanetary medium and their relation to the propagation and 
acceleration of cosmic rays. Some related work by others will be referred 
to, but no attempt is made to represent all such work. 

The discontinuities which are observed near 1 AU are hydromagnetic 

(MHD) discontinuities - hydromagnetic because the magnetic pressure is 

comparable to the thermal pressure near 1 AU. The thermal pressure has 

been uncertain until recently because of ignorance of the electron 

temperature T , but this situation has been remedied by both indirect and 

direct determinations of T , discussed in Section II. With this result, 

e 

the important fact that the solar wind is a P«1 plasma is now established. 

Three kinds of MHD discontinuities are important to cosmic ray physics: 
tangential discontinuities, rotational discontinuities, and shocks. 
Tangential and rotational discontinuities are important in theories of 
cosmic ray propagation. Shocks are of special interest because they can 
accelerate cosmic rays. 

Tangential and rotational discontinuities interact in distinctly 
different ways with cosmic rays. Despite their fundamental physical 
difference, however, it is difficult in practice to distinguish between a 
rotational and tangential discontinuity. This has led to a controversy 
concerning the basic questions of the relative abundance of these two 
types of discontinuity and their contribution to the power spectrum. The 
controversy seems to have been resolved, however. These results are 
presented in Section III together with the implications concerning cosmic 



ray propagation. 

A series of interesting experimental and theoretical papers has 
recently appeared concerning the acceleration of particles by forward 
fast shocks. New kinds of shocks have also been discovered. Some of 
this work is discussed in Section IV. 
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II. P at 1 AU 

One of the most important characteristics of the solar wind is the 
ratio of the thermal pressure to the magnetic pressure, 

nk(T +T ) 

P - 

B /8tt 


Measurements of the density n, proton temperature 



and magnetic field 


intensity B have been available for years, but the electron temperature 


T has remained obscure until recently. Now both indirect and several 
6 

direct measurements of T have been obtained. The indirect measurements 

e 

(Burlaga, 1968; Burlaga and Chao, 1971; Ogilvie and Ness, 1969) are based 
on the pressure balance condition at tangential discontinuities, 


8tt + nk <VV 


The direct measurements are from electrostatic analyzers of Bame, Ogilvie, 
and Serbu. All of these determinations give essentially the same result, 

T - (1.5 ± .5)x10 5 °K 

The electron temperature is nearly constant. It changes little with 

time. Unlike the proton temperature, T is independent of the bulk speed 

e 

(Burlaga and Ogilvie, 1970a) (see Figure 1) and of gradients in the bulk 
speed (Burlaga et al., 1971). Burlaga et al. (1971) suggested that the 
near constancy of T is due to the high thermal coniictivity of electrons 
near 1 AU this has been supported by more a detailed study of Hundhausen 
end Montgomery (1971). 

5 

Using T g = 1.4x10 °K, Burlaga and Ogilvie (1970b) found that the most 
probable value of (3 was unity in the period June-Movember, 1967 (Figure 2). 
Thus, the magnetic pressure and the thermal pressure tend to be equal, 
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and so p » l, For this reason, one must use magnetohydrodynamics rather 
than hydrodynamics to describe the solar wind. The form of hydromagnetic 
theory which is applicable to the solar wind near 1 AU has been presented 
by Burlaga (1971a). 








Mamm 


'll MliiTrilVBiiBrtMMiliiii 


II. Tangential and Rotational Discontinuities 

A. Scattering properties of RD's and TD f s 

Although several people have studied the interaction of cosmic rays 
with rotational and tangential discontinuities, most of these computations 
remain in their files; relatively little has been published. The theory 
is straightforward but the results are complex and the application to a 
realistic propagation situation hi even more complicated. Here we shall 
simply say a few words about the important differences between scattering 
by RD's and TD's. 

Consider the interaction of a charged particle with a tangential 
discontinuity. The discontinuity may be viewed as a surface in space. 
Generally, such surfaces are appreciably bent on a scale of .01 AU, at 
least near the earth, (Burlaga and Ness, 1969), but let us consider a 
plane surface. The magnetic field is parallel to the surface, but otherwise 
arbitrary. Consider the simplest case of oppositely directed fields on 
the two side?) of the surface (Figure 3a), If the field is uniform, a partible 
will move along the surface, perpendicular to JB. A tangential discontinuity 
can thus act as a highway in space, transporting particles perpendicular 
to JB at a rapid rate nearly equal to the particle speed. In the absence 
of such a discontinuity, the particles would be constrained to move parallel 
to JB. Generally the fields are not anti-parallel, but make some small angle 
to one another (Burlaga, 1969a). The same principle applies however; the 
result is that the drift rate is lower in this case and the trajectory is 
more complicated. 

In practice fluctuations, gradients, and the curvature of the surfaces 
will tend to remove particles from the surface; in addition one must consider 
how they arrive at a surface. This problem is full of complications. It 
is being studied by Fisk and Schatten who suggest that TD's near the sun 




may provide paths to carry particles rapidly from the east limb to the 
west limb. 

Although TD's can transport particles , they do not act as scattering 
centers. A region containing only TD's and uniform fields would not act 
as a diffusing medium. 

Now consider scattering by a rotational discontinuity. Such a 
discontinuity can be viewed as a kink in the magnetic field lines figure 3b). The 
problem was treated by Parker (1963). When particles move froip a weak 
field to a strong field, some are reflected; the transmission coefficient 
is f * cos 0 where 0 is the angle between B and the direction of the particle 
flux. When particles move from a strong field to a weak field some might 
be reflected, but they return and ultimately all are transmitted. 

The result is that some rotational discontinuities can act as efficient 
scattering centers for cosmic rays. Thus, a region containing numerous 
rotational discontinuities could act as a diffusing medium. 

Quenby (1971) has suggested that rotational discontinuities may have 
been the dominant scattering centers in a diffusive cosmic ray event that 
he analysed. He also suggested, based on indirect evidence, that TD's 
and RD's occur equally frequently, but the results in the next section do not 
support this. 

B . Relative number of TD's and RD's 

As discussed above, TD's and RD's Interact very differently with cosmic 
rays, so it is essential to determine the relative number of tangential 
and rotational discontinuities. This has been a controversial subject. 

Smith et al. (1970) stated that most interplanetary discontinuities are 
rotational, Davis (1970) and Quenby (1971) suggested that tangential and 
rotational discontinuities occur with nearly equal frequency, and Ness 
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meaningful discussion of this subject requires an operational definition 
of a discontinuity. We shall use the definition given by Burlaga (1969a) 
for a "directional discontinuity", which, with some simplification, is 
a change >30° in the magnetic field direction which occurs in <30 sec as 
a result of the motion of the discontinuity past the spacecraft. (Some 
directional discontinuities are shown in Figure 4). This definition includes 
the discontinuities discussed by Ness, it is nearly equivalent to that for 
the discontinuities analyzed by Siscoe et al. (1968), and it specifies the 
kinds of discontinuities referred to by Smith et. al. (1970). 

Recent results of Burlaga (1971) indicate that most directional 
discontinuities are tangential. The argument is as follows: If directional 

discontinuities are rotational, then the charge in the velocity across a 
discontinuity would be related to the change in the magnetic field. 
Specifically, Hudson (1970) showed that 



fhis can be written 



(3) 

(4) 


A depends on the thermal anisotropy of the plasma and is .9 ± .1 for 
typical solar wind conditions near 1 AU. 
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Applying the above test to the 200 diecontimiitiee which peased 
Pioneer 6 in the period Dec. 18-25, 1965, Burlaga (1971b) found the 
reeult in Figure 5. Clearly, (3) is not satisfied. Instead, the most 
probable value of AV, is aero. Thus, most directional discontinuities 
are not rotational discontinuities. Since they do not have signatures 
characteristic of shocks, we infer that most of the discontinuities are 
tangential. 

Directional discontinuities are not isolated features like shocks. 
They form a complicated network of surfaces which underlie; the basic 
structure of the solar wind, as suggested by the artistfs drawing in 
Figure 6, from Burlaga (1971c). They occur throughout the region 
between 8 AU and 1 AU (Burlaga 1971b) and are separated by « .01 AU 
(Burlaga, 1969). This separation is comparable to the mean free path 
obtained from models of solar cosmic ray events; since directional 
discontinuities are mostly tangential and thus do not scatter cosmic 
rays appreciably, this relation is probably coincidental. 

What then is the role of directional discontinuities in cosmic 
ray transport? Although they do not act as scattering centers, they 
must be considered in cosmic ray diffusion theories for the reason 
discussed in the next section. 

C. Contribution of Discontinuities to Power Spectra 

Sari and Ness (1969) showed that at certain times directional 
discontinuities dominate the power spectrum, i.e. Th® power spectrum 
calculated from the discontinuities alone is identical to the measured 
power spectrum. This presents a problem for cosmic ray propagation 
theories: calculating the diffusion coefficient from the power spectrum 
using the standard scattering theory would imply appreciable scattering 

. ... . ... . ... .. . ..... ^ .... *-*- - ■- ■ - ' •, || |j) ,i 
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and diffusion in this case, whereas in fact the particles would not 
diffuse appreciably since most of the power is due to tangential 
discontinuities. 

Clearly, one must determine how much power discontinuities con- 
tribute to the spectrum before he calculates diffusion coefficients 
from a theory which assumes scattering by waves. It is important, 
therefore, to determine how frequently the spectrum is dominated by 
discontinuities and what fraction of the power is due to discontinuity 
at other times. This problem is %eing studied by Sari at Goddard. 

His investigation is not complete, but the preliminary results indicat 
that discontinuities dominate the spectrum only occasionally, when the 
solar wind is quiet (low speed and temperature) . At other times, dis- 
continuities make a similar contribution to the power spectrum, but in 
addition waves, non-linear fluctuations, and various types of static 
structures contribute to the spectrum in still greater amounts. An 
example of the type of static structures chat are present is shown in 
Figure 7 from Burlaga and Ogilvie (1970b). Any analysis of cosmic 
ray propogation which uses power spectra should consider the kinds of 
features that are contributing to the spectrum during the interval of 
interest. This changes appreciably from day to day. 

IV. Shocks 

A. Existence of various types of shocks . Shocks can be classifie 
in several ways (e.g. ; see Burlaga, 1971a)* Here we shall consider 2 
types of shocks - fast shocks and slow shocks. Relative to the solar 
wind, each of these types can propagate either away from the sun 
(forward shocks) or toward the sun (reverse shocks). One can thus 

speak of four kinds of shocks: 
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forward fast shock 
reverse fast shock 
forward s i ow shock 
reverse slow shock 

The reverse shocks move away from the sun even though 
they move toward the sun relative to the plasma, because they are 

convected outward by the solar wind. An important difference be- 

tween a fast shock and a slow shock is that the magnetic field 
intensity increases across a fast shock but decreases across a 
slow shock. 

» 

The existence or fast forward shocks is well known. Currently, 

interest is centered on the shape of the shock surface?, the inter- 

action of shocks with other discontinuities, and the effects of the 
thermal anisotropy on the jump conditions. Lepping and Argentiero 
(1971) have developed a technique for accurately calculating shock 
normals which will be of use in studying the first two problems above. 
Lepping is currently studying the interaction of a shock with the 
bow shock. He has also shown that the anisotropy does not 
significantly change the jump conditions, i.e. the use of the Rankine 
Hugoniot conditions for an isotropic medium is a good approximation 
for fast shocks. 

The existence of a reverse fast shock in the solar wind was 

demonstrated by Burlaga (1970) using Explorer 34 plasma and magnetic 

field data obtained by Ogilvie and Ness, respectively. The basic 

observations are shown in Figure 8 together with a plot of the 
2 

pressure P = B / (8rr) + nk(T + T ). n, T, and B decrease with time 

© 
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because one sees the back side of the shock first. This was, in fact, a 
perpendicular shock, i.e. and Bg were perpendicular to the shock normal 

and parallel to one another. As required for a perpendicular shock, 

11 ^ 

B-/B p = — . The shock speed and direction were determined using simultaneous 
observations from Ness 1 magnetometers on Explorers 33, 34, and 35. The 
shock normal was found to be in the ecliptic plane and perpendicular to 
the spiral direction (~45° from the earth-sun line). It moved toward the 
sun, relative to the solar wind, at a speed V - 141 km/sec. It was verified 
that the Rankine-Hugoniot conditions were satisfied for this shock. The 
origin of this shock is a mystery. The shock appears ahead of a high speed 
stream, suggesting it might be due to a stream-stream interaction. This 
should produce both a forward and a reverse shock. However a forward shock 
was not seen. 

Evidence for forward slow shocks in the solar wind was first published 
by Chao and Olbert (1970). Additional evidence for such a shock was 
presented by Burlaga and Chao (1971) (see Figure 9). The latter shock 
occurred at ~1423 UT on Jan. 20, 1966. Note that the shock does not appear 
as a sharp discontinuity in the magnetic field data, and that there are 
large fluctuations in B near the shock. These characteristics t §ommon 
to all the slow shocks that have been observed, make the identification 

jf slow shocks rather difficult. The procedure used by Burlaga and Chao (1970) 
;as to find a solution to the Rankine-Hugoniot equations which was consistent 
rith the observations within the uncertainties due to the measurements and 
luctuations, and which satisfied other necessary conditions for a slow 
hock. The solutions to the R-H equations are shown by the horizontal lines 
n Figure 9. Clearly they are consistent with the observations, but one 
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would have hoped for smaller fluctuations in the data. Additional evidence 
supporting the identification of the discontinuity in Figure 9 as a slow 
shock is the following: (a) predicted values of Vg are in good agreement 
with the observed values, (b) the flow speed normal to the shock did 
decrease across the shock (from 37 km/sec to 27 km/sec) , (c) the Alfven 
Mach number was less than 1 on both sides of the shock (.9 ahead and .8 
behind), and (d) the slow mode Mach number was greater than 1 ahead of the 
shock (1.3) and less than 1 behind it (.8), as required by the theory of 
slow shocks. 

A reverse slow shock was also found in the Pioneer 6 data by Burlaga 
and Chao (1971). The observations and a solution to the Rankine-Hugoniot 
equations are shown in Figure 10. Here too it was found that (a) predicted 
values of Vg are in agreement with observed values, (b) the flow speed 
decreased across the shock (from 29 km/sec to 23 km/sec), (c) the Alfven 
number was <1 on both sides of the shock (.9 ahead, .8 behind), and (d) the 
slow inode Mach number was >1 ahead of the shock (1.2) and <1 behind it (.8), 

Thus, there is evidence for all 4 of the types of shocks listed at 
the beginning of this section. Ivanov (1970) predicted still another kind 
of shock, which he called a "rotational discontinuity" because it should 
propagate at nearly the Alfven speed. The coplanarity theorem applies 
for such a shock, but unlike most shocks, n^^g. The magnitude of B does 
change, and the change in entropy is due to a change in the anisotropy. 
Ivanov suggested that most of the discontinuities in Burlaga (1969b) were 
such shocks, but a more detailed study by Burlaga (1971c) shows that this 
is not the case. Burlaga and Chao (1971) show that discontinuities with 
B^Bg and n^=ng seldom if ever occur. 




' it i 
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Given the existence of fast and slow shocks, one may now ask how 
they relate to cosmic rays. In the case of fast forward shocks there is 
a very fundamental relation which is discussed in the next section - 
shocks can accelerate cosmic rays. Whether or not reverse fast shocks and 
slow shocks can accelerate particles remains to be determined. This 
problem has not yet been investigated, but it promises some interesting 
results . 


B. Proton acceleration by forward fast shocks. 

Ogilvie and Arens (1971) showed that increases in the flux of 1-10 Mev 
protons are sometimes associated with forward fast shocks (Figure 11). 

Similar increases were observed by Armstrong et al., (1970) who found that 
heavy particles (Z>2) are accelerated as well as protons but electrons 
are not accelerated. Such increases occur only when there is an appreciable 
flux of 1-10 Mev particles present before the arrival of the shock. The 
peak flux is «10 times this ambient flux. The spectra of the ambient flux tends 
to be rather steep. These increases are brief («10 min), and are thus 
distinct from the kind of shock-associated events discussed by Vernov et 
al. (1970). 


It is generally agreed that such increases are due to the acceleration 
of ambient energetic particles. Acceleration lifts the more abundant 
low energy particles above the detector threshhold and thus gives an 
apparent increase in the flux of particles at energies near the threshhold. 
The flux increase should be greater for steeper spectra, for a given 
acceleration mechanism; the observations of Ogilvie and Arens show this 
effect. 

The nature of the accelerating mechanism is controversial. Axford and 
Reid (1963) had suggested that particles would be accelerated if they were 




• 






trapped between the earth* s bow shock and an interplanetary shock moving 
toward the earth. This mechanism requires only that there is a certain 
probablility for reflection of particles by the shock and that the magnetic 
field lines which guide the particles should intersect both shocks. The 
latter condition was met for all except perhaps one of the five shock 
associated increases observed by Ogilvie and Arens, so their results were 
interpreted as support for the Axford-Reid model. They did not exclude 
other accelerating mechanisms, how ever. 

An alternative accelerating mechanism was proposed by Fisk (1971). 

His idea is that particles are retained near the shock by diffusing in 
the magnetic field fluctuations which are ahead of the shock. In this 
model, the fluctuations play the role of the bow shock in the Axford-Reid 
model . 

The flux of reflecting low energy particles at the shock is 


S 

n 


V U - I V != (TU ) 
s r 3 s ST r' 


(5) 


where V is the shock speed, T the particle kinetic energy and U is the 
s r 

differential number density of reflecting particles which is proportional 
to the upstream density of particles, U. The density of particles U is 
determined by solving the convection diffusion equation 


SU 
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where V’ and K* are the wind speed and diffusion coefficient ahead of the 
shock. A similar equation describes particles behind the shock, U*'. The 
solution for large times is 
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where T| ~ r-V B t (V being the shock speed), V 0 « V -V' , and 0(l+2p,)/3 

(|i being the spectral index, U * A(T/T It is assumed that the 

shock is nearly a perfect reflector for low energy particles* Equation (7) 

has several free parameters: K' , T 1 , p>, V and V'. Only 2 of these 

parameters can be determined by fitting the data. Fisk uses typically 

measured values for V g and V' , and he then finds pairs of value of \i and 

which give the measured peak intensity at the shock. K' is 

determined uniquely. An example of such a fit for the Nov. 29, 1967 

18 2 

event is shown in Figure 12. It gives K' ~ 10 cm /sec, or K (parallel 
19 2 

to B) ~10 cm /sec. Clearly, Fisk’s model can fit the observations and 

gives reasonable values of K' , but the fit is not unique. 

If Fisk’s model is correct, one expects to find shock-associated 

increases far from the earth. On the other hand, if the Axford-Reid 

mechanism is the only one that operates, one should not find shock- 

associated increases far from the earth. Well, such an increase has been 

8 

observed by a deep space probe, Pioneer 8, 1.5x10 km (^1 AU) from the earth 
by Palmeira et al. (1971). Strong support for Fisk’s mechanism! 

But further complications were revealed in a paper by Singer (1970). 

For example, he observed strong anisotropies in the energetic particles 
near the shock, which seems inconsistent with Fisk's idea of diffusion. 

He also found that the intensity maximum can occur 5-10 min after the 
shock. Singer believes that the gain of a factor of ~5 occurs in a single 
encounter rather than by multiple encounters as in the Fisk and Axford- 
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Reid models. 
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Clearly, much remains to be learned about acceleration by shocks. 

V. Summary 

This brief review has considered several distinct, but closely related 
topics: the electron temperature and value of 3 in the solar wind, the 
relative number of tangential and rotational discontinuities and their 
relation to the theory of cosmic ray propagation, the existence of various 
types of shocks, and the acceleration of particles by fast forward shocks. 
Rapid progress in these areas has been made in recent years due to the 
availability of good data and to the development and application of 
relevant theories. Several new problems have also come into focus: What 

is the thermal conductivity of electrons? What are the origins of the 
various discontinuities? How does one describe cosmic ray propagation in 
a medium where the power is due to structures other than linear waves? 

How are particles reflected by a shock? Which of the several accelerating 
mechanisms is dominant in various circumstances? Hopefully, the answers 
to these questions will appear within the next several years. 
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Temperature speed relation. Explorer 34 observations show 

a linear relation between v/Y" and U, where T is the 

v p p 

proton tempera ture and U is the bulk speed. Observations 

from other experiments fall on the same line. The electron 

temperature observations, although still sketchy, suggest 

that T is independent of U; this is indicated by the 
e 

horizontal line. Some theoretical models are also shown. 

See Burlaga et al (1970) for the details. 

Distribution of (3 based on Explorer 34 data. The most probable 
value is 1.0 sis 1. 


Tangential and rotational discontinuities. The top illustration shows 
how a tangential discontinuity (viewed edge on) transports a 
cosmic ray along it surface when and are antiparallel. 

The lower illustration shows a special kind of rotational dis- 
continuity. 

Some directional discontinuities. 

Distribution of AV/Q for % 200 directional discontinuities. 

Since the peak does not occur near .9, most of these discontin- 
uities are not rotational. 

Simplified view of 3 discontinuity surfaces and magnetic fields 
between them, illustrating how a .05 AU segment of the solar 
wind might look. 
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Fig, 7 The magnetic pressure P^ and thermal pressure P^ tend to 

be anticorrelated on a scale of .01 AU tending to keep the 
total pressure P^, constant on that scale, even though P^ 
changes on a larger scale. The anticorrelation suggests sta- 
tic features. These will contribute to the power spectrum 
of B. 

Fig, 8 A reverse fast shock. 

Fig. 9 A forward slow shock. 

Fig. 10 A reverse slow shock 

Fig. 11 This shows an increase in the intensity of ~ 1 MeV cosmic 
rays ahead of a shock. 

Fig. 12 This shows Fisk f s theoretical fit to the observations shown in 
Figure 11. 
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